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The absence of locahzation of impurity-induced low-energy quasiparticle states in a two-dimensional 
d-wave superconductor is argued for any amount of disorder in the limit of unitary scatterers. This 
surprising result follows from the fact that a unitary impurity produces a marginally-bound state at 
zero energy which decays as a power-law along the nodes of the d-wave energy gap. Consequently, for 
finite density of impurities, the impurity-induced states are coupled by long-range overlaps yielding 
extended quasiparticle states below a characteristic energy scale tOc- Simple scaling arguments 
suggest that cjo oc e~'^^"'"p, where riimp is the impurity density and c is a positive constant. 
PACS numbers: 74.25. Jb, 71.27.+a 



When impurity scattering violates the symmetry of the 
superconducting condensate, the superconducting energy 
gap is depleted and idtimatcly destroyed by impurities. 
This happens, for example, in 5-wave superconductors 
with magnetic impurities [1,2]. Some of the recent exper- 
iments support d-wave symmetry as a pairing channel in 
the cuprates [3-5], predicted by the spin-fluctuation the- 
ory [6-8] as well as by the Hubbard model [9]. For this 
state, even scalar impurities act as pair-breakers, produc- 
ing a finite lifetime for quasiparticles near the nodes in 
the gap, and a finite density of states at low energies [10]. 
As a result, the measured low temperature properties of, 
for example, YBa2Cu3 07 [3], Bi2Sr2CaCu2 08 [11], and 
Lai,86Sro,i4Cu04 [12] display a remarkable sensitivity to 
the presence of impurities. 

Recently, the role of imperfections in d-wave super- 
conductors has been considered by Lee [13], Hirschfeld 
and Goldenfeld [14], and subsequently by others [8,15,16]. 
Qualitatively these low-energy quasiparticle states be- 
have as a disordered Fermi liquid with a renormalized 
density of states. Impurity scattering in a disordered 
s-wave superconductor can lead to localization of low- 
energy states, as was first pointed out by Ma and Lee [17] 
and Mackawa and Fukuyama [18]. The same effect should 
also, in principle, take place in disordered d-wave super- 
conductors [13,16]. Lee showed that in two-dimensional 
(2D) disordered d-wave superconductors these states are 
subject to weak-localization corrections and thus are lo- 
calized. Localization of quasiparticles was advocated as a 
way to recover some .s-wave features in a d-wave system. 
Here, we argue that the conclusion regarding localiza- 
tion found in Refs. [13,16] is based on a simplified picture 
which does not take into account the long-range nature 
of hopping between impurity states. It is this long-range 
hopping which makes our approach different and which 
is lost upon trivial averaging over impurity positions. 

In obtaining the above results, impurity-averaged 
quantities were expressed in terms of ensemble-averaged 
single-particle Green's functions which were calculated 
within the coherent-potential approximation (CPA) [19]. 
Note that, in the CPA, the quasiparticle self-energy 



S(ci;) = ninipT{io) is approximated by a single impu- 
rity T-matrix; nimp is the impurity density. Spatial 
anisotropy of impurity states, which is important for a 
d-wave superconductor, is lost upon averaging. More- 
over, the CPA is not adequate for studying fine details of 
impurity bands [19]. 

In this Letter, we consider the problem of localization 
of quasiparticle states at low energies, starting from a 
single-impurity solution as a basis [20]. Our main result 
is the conclusion that the weak-localization theory does 
not apply to localization of impurity-induced quasiparti- 
cle states in 2D d-wave superconductors because of the 
long-range interactions between these states [21]. We ar- 
gue that, in order to address this question, one has to 
implement an approach where the single-impurity prob- 
lem is solved first. This approach was used by Shiba and 
Yu to study the eff'ect of magnetic impurities in an s-wave 
superconductor [22,23]. They showed that a quasiparti- 
cle bound state is formed in the energy gap as a result 
of multiple scattering. A similar consideration of scalar 
impurities in a p-wave superconductor was done by Bucli- 
holtz and Zwicknagl and by Stamp [24]. These mid-gap 
states eventually form an impurity band in conventional 
superconductors. Notably, scalar unitary impurities cre- 
ate largely analogous impurity states in d-wave supercon- 
ductors [20]. 

Below, we first argue that, due to the long-range over- 
laps between impurity states, disorder scattering does not 
localize quasiparticle states at w = and in its vicinity. 
Along the directions of the vanishing energy gap (the di- 
agonals of the square lattice for a d^2_y2 gap function), 
the impurity wave functions decay as 1/r. This in turn 
leads to the 1/r power-law overlaps between impurity 
states and to a novel network of strongly coupled im- 
purities, yielding extended impurity states (see Fig. 1). 
From the point of view of the localization theory, our 
impurity problem belongs to a new class where for any 
amount of disorder long-range hopping delocalizes im- 
purity states. Second, we speculate that the nature of 
low-energy states in a disordered d-wave superconductor 
is qualitatively different from the simple extended states 
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FIG. 1. (a) Three impurities are shown at 1,2, and 3. For 
the impurities at sites 1 and 2, the long-range 1/r tails of the 
impurity states along the nodal directions of the energy gap 
A(</i) are depicted by two parallel lines. Strong overlap be- 
tween these states leads to a resonance behavior, if the on-site 
energies ei and £2 are close enough. The exponentially small 
overlap between states at sites 1 and 3 will not produce a res- 
onance at small impurity densities, ^o^imp ^ 1, and are omit- 
ted reducing the original problem to a locaUzation problem 
on a network of strongly coupled impurity states. The ring 
(R, R + clR) denotes the area which specifies the probability 
PresiR^dR) [see Eq. (4)] for having another impurity resonant 
with the impurity at site 1. (b) The density of states N of 
a disordered d-wave superconductor. Crossover from the lin- 
ear dependence N{lu) oc o; to the impurity-dominated regime 
occurs at a;* oc n^^^ for unitary scatterers [13,14]. Due to 
the strongly-interacting impurities, the low-energy quasipar- 
ticle states are delocalized. At higher energies, the impurity 
states should form a band in which weak localization should 
prevail [13] and there should be an "inverse" mobility edge at 
ujc (oc e-'=/"™p). 

predicted in any theory with averaging over random im- 
purity ensembles. 

Originally, the effect of long-range hopping on lo- 
calization was considered by Anderson [26] for ran- 
domly distributed impurities in D dimensions with 
the Vij ~ hopping interaction (c/., RKKY- 

interacting spins). He showed that, for e > 0, there is 
a localization transition; for e < 0, states are always de- 
localized; and e = represents the marginal case with 
power-law decaying states. Subsequently, an analogous 
problem of localization of phonon modes in disordered 
solids was considered by Levitov [27]. We will use a sim- 
ilar method, modified to take into account strong spatial 
anisotropy of impurity states. 

First, consider overlaps between two impurity sites. 
For a unitary-scattering impurity, the bound state is well- 
defined and is located at zero energy. The crucial point in 
our approach is that, for a 2D d-wave superconductor, the 
impurity bound-state wave function is highly anisotropic 
and forms a cross-shape state with four tails along the 
diagonals of the square lattice [20]: 



^imp(T', <?!') oc sinkpi 



where (,{(!>) = hvF/\A{(f>)\ = ^o/|cos2(?i| and n = 0, 1,2,3 
correspond to the four nodes of the gap, along which the 
tails are slowly decaying. We will consider, having in 
mind strong planar anisotropy of cuprates, a 2D rf-wavc 
superconductor: A{(f)) = Aocos2(/i, where Aq is the gap 
amphtude and (j) is the planar angle. 

In the tight-binding approximation, which is valid at 
low impurity densities, ^oWimp <C 1, the hopping matrix 
element between two impurity sites and rj is given 
by Vjj = G(°^(ri — rj,0), where &°\t,uj) is the unper- 
turbed Green's function in the superconducting state (for 
details, see Ref. [28]). We shall work in Nambu space 
where "pj is the Nambu spinor describing a quasiparti- 
cle(hole) created in the ith impurity state, Eq. (1), and 
Vij is a matrix in this space spanned by the Pauh matri- 
ces, Ta {a = 1,2,3) [29]. The Hamiltonian is: 
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where the random distribution of on-site energies €{ is 
assumed to be uniform: P(ej) = l/W, for < < W, 
and zero otherwise. Since we are working with Nambu 
spinors, the on-site energy is a matrix too. Strictly speak- 
ing the impurity states are well defined only at cj = 0; by 
assigning the on-site energies a distribution, these states 
will be hybridized with the continuum of quasiparticle 
states [for a clean d-wave superconductor, N{(jj) oc (jj]. 
We shall ignore this small hybridization and consider only 
the direct overlap of the impurity-induced states. 

Using Eq. (1), one can show that the overlap decays 
exponentially (these sites will be called weakly overlap- 
ping), unless two impurities are connected by a radius 
vector r with the angle (p = 7r/4 -I- n7r/2, in which case 
the tails of the impurity states are strongly overlapping 
[e.g., the impurities 1 and 2 in Fig. 1(a)]. These sites 
are called strongly overlapping. In this case, the hopping 
matrix elements decay as 1/r: 
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where kp is the Fermi momentum, Vq is the dimension- 
less strength of the matrix element, and 6ij is a phase 
shift, which we will ignore hereafter. Both the fi and 
fa components of the overlap integral are present in gen- 
eral. However, the fi component of Vij is zero along the 
diagonals [28]. 

To establish delocalization, we show that for a given 
impurity site there exists a (large) distance R within 
which this impurity will always find another strongly 
overlapping impurity. The long-range hopping will pro- 
vide a strong resonance between these sites and the wave 
function will be delocalized [26,27]. This should be con- 
trasted with the weak-localization theory which predicts 
that all quasiparticle states are locahzed in 2D. 

We calculate the probabihty density for a given impu- 
rity state at site i with energy ej to have a resonance with 
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another impurity state located at a distance R, as shown 
in Fig. 1(a). For any two sites with on-site energies 
and ej, the resonance will occur only for a large enough 
hopping, \Vij \ > — fjl- The fs component of F, 1'3(-R), 
decays with the distance, whereas the phase space for 
possible resonances is growing with the distance and the 
balance for the resonance probability will depend on the 
power of the decay of V^iR) [26,27]. For l/i?-decay with 
essentially one-dimensional motion along the tails of the 
wave function, we find: 

P,MdR) = a^ln{R/^o)\smkFR\, (4) 

where a = i^QTiimpVohvp /{W^o) is a small dimension- 
less parameter. This small parameter allows us to ig- 
nore higher order simultaneous resonances. The average 
number of resonances / Pi-es{R,dR) = N^es, is divergent 
with distance and, therefore, we arc led to consider the 
problem of strongly overlapping random impurity states, 
where the tails of the impurity wave functions give rise 
to the network of strongly coupled impurity states [30]. 

To derive Eq. (4), we omit weakly overlapping sites 
because, for low densities ^oKimp ^ 1, their contribution 
is exponentially small. The inclusion of the weakly over- 
lapping sites will only help delocahzation we are set to 
prove. Next, we define the angular size of the tails at 
large distances, as seen from the impurity site. Matching 
two asymptotics of the impurity wave function, Eq. (1), 
one finds that the exponential asymptotic matches the 
power law at the angle = 7r/4 ± S{R), where 

6{R) = ^HR/^o). (5) 

The logarithmic factor indicates that the tails are con- 
tinuously broadening at large distances. The probabil- 
ity density is given as Pyes{R, dR) = J d4)Pres{R, dR, d(j)) 
with angular dependent probability being nonzero only 
for \<1) - 7r/4: + n7r/2\ < 2S{R) and zero elsewhere. This 
constraint results from the strong spatial anisotropy of 
the overlaps, as shown in Fig. 1(a). We find that 

\V-i(R)\ 

Pres {R,dR,d(l>)= n,^pRdRdcj>^-^^ . (6) 

Combining Eqs. (5) and (6), we finally obtain Eq. (4). 

Consider the average number of resonances for a 
given site inside the circle of radius R: N^^siR) = 
J PTes{P,dR) (X a\n^ R/^o- Diverging number of reso- 
nances indicates delocahzation. Another interesting re- 
sult is the In^ R/^o behaviour of the probability den- 
sity. 1/R decay of the hopping elements along the one- 
dimensional tails would lead to the Ini? growth of N^es- 
However, the second logarithmic factor comes into play 
due to the cut-off angle S{R) and the problem is "su- 
percritical" instead of being simply marginal. To illus- 
trate this point, we calculate the probability of having 



no resonances -Pnores(-R) between a given site and any 
other sites inside a circle of radius R. By drawing con- 
centric circles Ri {i = 1,...,N), Ri+i — Ri = p with 
Co p «C i? around the chosen impurity and calculating 
-Pnores(-R) = HJl " -^res (-Rj , p)] , wc find log-uormal decay 
of no-resonance probability instead of the expected slow 
power-law: 

Pnores 

(J?) ocexp[-0(l)aln'(J?/Co)]. (7) 

One can estimate the average minimum distance between 
resonating sites as Rc (x J RdR exTp[—0{l)a\n^ R] oc 
exp[0(l)/Q!], which generates an energy splitting coc = 
hvp/Rc. We interpret uic as the characteristic energy 
scale of the delocalized states at < uJc', see Fig. 1(b). 

There is a caveat to the interpretation of the above cal- 
culation. The probability of pairwise resonances, Eq. (4), 
does not take into account the effect of closed loops, 
which provide hacks cattering, responsible for the local- 
ization. We believe, however, that loops are irrelevant at 
small a. It is easy to show that a closed path should con- 
tain at least four sites, since at each impurity site particle 
can only turn by an angle ±n7r/2 in order to move along 
the tails. Due to this geometric constraint the probabil- 
ity of quartic resonances is proportional to and does 
not diverge at large R [28]. Therefore, at large distances 
pairwise resonances dominate. Based on this argument, 
we may conclude that inclusion of the loops should not 
change our basic result. In the opposite hmit of no on- 
site disorder W ^ Q {a ^ \), the problem is in the 
strong coupling limit and multiple resonances are impor- 
tant. In this limit, the pairwise-resonance approximation 
is invalid and another approach is necessary. Nonethe- 
less, it is reasonable to assume that the impurity states 
are delocalized again in this limit due to the strong over- 
laps between impurity states. However, oJc becomes of 
the order of or larger and our approach breaks down 
because the hybridization with the continuum is a new 
relevant feature in the problem. 

Finally, we comment on the nature of the delocalized 
states. In principle we are facing two possibilities: («) the 
states are simply extended with nondecaying probability 
density in the whole sample. In the view of the highly 
nontrivial form of Pres(^), this seems unhkely, although 
it cannot be excluded at a moment, {ii) The extended 
impurity states are "critically" locahzed, e.g., the impu- 
rity wave function envelope decays as l/r*^ with some 
index < < 1. This could happen if the "dressing" of 
single impurity wave functions due to resonances is not 
sufficient to delocalize the states completely. At present 
we cannot distinguish between these two possibilities [31]. 

The results presented here are sensitive to the spe- 
cific form of the energy gap close to the impurity. While 
the energy gap may change near the impurity, symme- 
try considerations [32] suggest that the node structure 
of the energy gap is not affected by the potential scat- 
terers. However, purporting a more general situation. 
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suppose that a small local imaginary 5-wave component 
of magnitude is generated in the neighborhood of the 
impurity. Such a component will most likely cut off the 
power-law tails of the wave function 0imp(i"), causing 
it to decay as r~'^e~^/^ along the diagonals with large 
£ = hvp/^s ^ ^0- In this case, the impurity states will 
be locaHzed for f^riimp <C 1. Nonetheless, there should 

— 1/2 

exist an intermediate regime, <C n-^^^p ^ where the 
strongly overlapping impurity states along the diagonals 
would lead to delocalization. 

In conclusion, we have studied formation of impurity 
bands in d-wave superconductors. We argue that long- 
range impurity-impurity interactions modify the usual 
weak-locahzation results and lead to extended quasipar- 
ticle states below a characteristic energy scale lOc oc 
g-c/ni„p ^ positive constant). This result suggests 
that unitary scatterers do not cause activated behavior 
of quasiparticles in quasi-2D d-wave superconductors at 
low temperatures. 
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